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Abstract. We solve the time-dependent Dirac equation for a muon which is 
initially bound to a fissioning actinide nucleus. The computations are carried 
out on a 3-D cartesian lattice utilizing the Basis-Spline collocation method. The 
muon dynamics is sensitive to the nuclear energy dissipation between the outer 
fission barrier and the scission point. From a comparison with experimental 
data we find a dissipated energy of about 10 MeV and a fission time delay due 
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1. Introduction 

There are two different mechanisms that contribute to nuclear energy dissipation, 
i.e. the irreversible transfer of energy from collective into intrinsic single-particle 
motion: two-body collisions and "one-body friction" . The latter is caused by the 
moving walls of the self-consistent nuclear mean field. The role played by these 
two dissipation mechanisms in fission and heavy-ion reactions is not yet completely 
understood. In a pioneering article that appeared in 1976 Davies, Sierk and Nix |l| 
calculated the effect of viscosity on the dynamics of fission. Assuming that friction 
is caused by two-body collisions they extracted a viscosity coefficient /i = 0.015 
Tera Poise from a comparison of theoretical and experimental values for the kinetic 
energies of fission fragments. The corresponding time delay for the nuclear motion 
from the saddle to the scission point was found to be of order At = 1 x 10 -21 s. 
However, in one-body dissipation models the time delay is an order of magnitude 
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larger. Several experimental techniques are sensitive to the energy dissipation in 
nuclear fission. At high excitation energy, the multiplicity of pre-scission neutrons 
H or photons || depends on the dissipation strength. At low excitation energy, 
the process of prompt muon-induced fission jj]] provides a suitable "clock" . This 
process will be discussed here. 

After muons have been captured into high-lying single particle states they form 
an excited muonic atom. Inner shell transitions may proceed without photon emis- 
sion by inverse internal conversion, i.e. the muonic excitation energy is transferred 
to the nucleus. In actinides, the 2p — > Is and the 3c? — > Is muonic transitions 
result in excitation of the nuclear giant dipole and giant quadrupole resonance, re- 
spectively, which act as doorway states for fission. The nuclear excitation energy 
is typically between 6.5 and 10 MeV. Most importantly, the muon is still available 
following these atomic transitions (in the ground state of the muonic atom) and 
can be utilized to probe the fission dynamics. Eventually though, the muon will 
disappear as a result of the weak interaction (nuclear capture by one of the fission 
fragments). However, the nuclear capture occurs on a time scale of order 1CP 7 s 
which is many orders of magnitude larger than the time scale of fission. 

The prompt muon-induced fission process is most easily understood via a "cor- 
relation diagram" , i.e. one plots the single-particle energies of the transient muonic 
molecule as a function of the internuclear distance If there is a large amount 
of friction during the motion from the outer fission barrier to the scission point the 
muon will remain in the lowest molecular energy level lscx and emerge in the Is 
bound state of the heavy fission fragment. If, on the other hand, friction is small 
and hence the nuclear collective motion is relatively fast there is a nonvanishing 
probability that the muon may be transferred to higher- lying molecular orbitals, 
e.g. the 2pa level, from where it will end up attached to the light fission fragment. 
Therefore, theoretical studies of the muon-attachment probability to the light fis- 
sion fragment, Pl, in combination with experimental data can be utilized to analyze 
the dynamics of fission, and nuclear energy dissipation in particular. 

2. Theoretical Developments 

Because the nuclear excitation energy in muon-induced fission exceeds the fission 
barrier height it is justified to treat the fission dynamics classically (no barrier 
tunneling) . For simplicity, we describe the fission path by one collective coordinate 
R; the classical collective nuclear energy has the form 

E nuc = \B{R)R 2 + VUR) + E M (R). (1) 

We utilize a coordinate dependent mass parameter p] and an empirical double- 
humped fission potential Vr s (R) pl which is smoothly joined with the Coulomb 
potential of the fission fragments at large R. The last term in Eq. (Q) denotes the 
instantaneous muonic binding energy which depends on the fission coordinate; this 
term will be defined later. 
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To account for the nuclear energy dissipation between the outer fission barrier 
and the scission point, we introduce a friction force which depends linearly on the 
velocity. In this case, the dissipation function D is a simple quadratic form in the 
velocity 

E nuc (t) = -2D = -fR 2 (t). (2) 

The adjustable friction parameter / determines the dissipated energy; it is the only 
unknown quantity in the theory. 

For the dynamical description of the muonic wavefunction during prompt fis- 
sion, the electromagnetic coupling between muon and nucleus (— ej^A^) is domi- 
nant; the weak interaction is negligible. Because of the nonrelativistic motion of 
the fission fragments the electromagnetic interaction is dominated by the Coulomb 
interaction 

A°( M )=|dV^^. (3) 

The muonic binding energy in the ground state of an actinide muonic atom amounts 
to 12 percent of the muonic rest mass; hence nonrelativistic calculations, while qual- 
itatively correct, are limited in accuracy. Several theory groups have demonstrated 
the feasibility of such calculations (|,|7],[8| which are based on the time-dependent 
Schrodinger equation 

[~ V 2 - eA°(r, t)} ^(r, t) = ih^ </-(r, t). (4) 

Recently, we have developed a numerical algorithm to solve the relativistic problem 
on a three-dimensional Cartesian mesh The time-dependent Dirac equation 

for the muonic spinor wave function in the Coulomb field of the fissioning nucleus 
has the form 

flb(i) Tp(r,t)=iH^ VM), (5) 
where the Dirac Hamiltonian is given by 

H B {t) = -ihca-V + f3mc 2 - eA°(r,t). (6) 

Our main task is the solution of the Dirac equation for the muon in the presence of 
a time-dependent external Coulomb field A (r, t) which is generated by the fission 
fragments in motion. Note the coupling between the fission dynamics, Eq. (|l|), and 
the muon dynamics, Eq. (Js[) , via the instantaneous muonic binding energy 

Ep{R(t)) = y>(T,t)\Hn(t)\il>(T,t)) (7) 



which depends on the fission coordinate; the presence of this term increases the 
effective fission barrier height. 
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3. Lattice Representation: Basis-Spline Expansion 

For the numerical solution of the time-dependent Dirac equation (||) it is convenient 
to introduce dimensionless space and time coordinates 

x = r/A c A c = Ti/(m^c) — 1.87/m 

r = t/r c t c = A c /c = 6.23 x l(T 24 s (8) 

where A c denotes the reduced Compton wavelength of the muon and r c the reduced 
Compton time. For the lattice representation of the Dirac Hamiltonian and spinor 
wave functions we introduce a 3-dimcnsional rectangular box with a uniform lattice 
spacing Ax. The lattice points are labeled {x a ,yp, z 7 ). 

Our numerical algorithm is the Basis-Spline collocation method ju^. Basis- 
Spline functions Bf 1 (x) are piecewise-continuous polynomials of order (M — 1). 
These may be thought of as generalizations of the well-known "finite elements" 
which are a B-Splines with M — 2. To illustrate the method let us consider a wave 
function which depends on one space coordinate x; we represent the wave function 
on a finite spatial interval as a linear superposition of B-Spline functions 

JV 

i>(x a )=Y,B?*(x a )c i . (9) 

i=l 

In the Basis-Spline collocation method, local operators such as the EM potential A 
in Eq. (|^) become diagonal matrices of their values at the grid points (collocation 
points), i.e. V(x) — ► V a = V(x a ). The matrix representation of derivative opera- 
tors is more involved |l(J. For example, the first-derivative operator of the Dirac 
equation has the following matrix representation on the lattice 

TV 
i=l 

where B' ai — [dBf 1 (x) / dx]\ x=Xa . Furthermore, we use the shorthand notation 
Bfji = B^ixp) for the B-spline function evaluated at the collocation point xp, 
and the inverse of this matrix is denoted by B 1 ^ = [i? -1 ]^. Because of the presence 
of this inverse, the operator D@ will have a nonsparse matrix representation. In the 
present calculations we employ B-Splines of order M = 5. Eq. (^) can readily be 
generalized to three space dimensions; in this case the four Dirac spinor components 
ip( p \p = (1, • ■ •, 4) are expanded in terms of a product of Basis-Spline functions 

^(x a ,y ,z^,t) = ]T Bf {x a )Bf (^)Sf (^)cg(t), (11) 

i.e. the lattice representation of the spinor wave function is a vector with N — 
4 x N x x N y x N z complex components. Hence, it is impossible to store -Hd m 
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memory because this would require the storage of N 2 complex double-precision 
numbers. We must therefore resort to iterative methods for the solution of the 
matrix equation which do not require the storage of Hp. 

We solve the time-dependent Dirac equation in two steps: first, we solve the 
static Coulomb problem at time t = 0, i.e. the muon bound to an actinide nucleus. 
This problem is solved by the damped relaxation method The second part of 
our numerical procedure is the solution of the time-dependent Dirac equation (^) 
by a Taylor-expansion of the propagator for an infinitesimal time step At. Details 
may be found in ref. H. 



4. Discussion of Numerical Results 

In the following we present results for prompt fission of 2 gjNp induced by the 3d — » 
Is muonic transition (9.5MeV). All results reported here are for a 3-D Cartesian 
lattice of size L x = L y = 67 fm and L z = 146 fm with N x X N y X N z = 25 X 25 X 53 
lattice points with a uniform lattice spacing Ax = 1.5A C = 2.8 fm. Depending on 
the value of the friction coefficient, we utilize between 1, 200 — 1, 900 time steps with 
a step size At = 1.5r G = 9.3 x 10~ 24 s. Typical production runs take about 11 hours 
of CPU time on a CRAY supercomputer or about 54 hours on an IBM RS/6000 
workstation. 

Fig. [l] shows the time-development of the muon position probability density 
during fission at a fragment mass asymmetry £ = Ah/Al — 1.10. As expected, 
the muon sticks predominantly to the heavy fragment, but for this small mass 
asymmetry the muon attachment probability to the light fission fragment, Pl, is 
rather large (20 percent). 

One might ask whether the muon will always remain bound during fission; what 
is the probability for ionization? To investigate this question we have plotted the 
muon position probability density on a logarithmic scale. 

In coordinate space, any appreciable muon ionization would show up as a "prob- 
ability cloud" that is separating from the fission fragments and moving towards the 
boundaries of the lattice. Fig. |2| shows no evidence for such an event in our numer- 
ical calculations. Hence, we conclude that the probability for muon ionization Pi on 
is substantially smaller than the muon attachment probability to the light fission 
fragment which is always clearly visible in our logarithmic plots, even at large mass 
asymmetry. From this we estimate that P; on < 10~ 4 . 

Fig. ^ shows that Pl depends strongly on the fission fragment mass asymmetry. 
This is easily understood: for equal fragments we obviously obtain Pl = 0.5, and for 
large mass asymmetry it is energetically favorable for the muon to be bound to the 
heavy fragment, hence Pl will be small. In Fig. |^ we also examine the dependence 
of Pl on the dissipated nuclear energy, Bdissj during fission. In our model, friction 
takes place between the outer fission barrier and the scission point. When the 
dissipated energy is computed from equation (0) we find an almost linear dependence 
of the muon attachment probability on Ed iss ; unfortunately, this dependence is 
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Fig. 1. Prompt muon- induced fission of 2 gjNp for a fission fragment mass asym- 
metry £ = Ah/Al — 1.10 at E* — 9.5 MeV. Shown is the muon position probability 
density at four different times during fission: t = 0, 6.5 x 1CP 21 s, 8.4 x 1CP 21 s, 
1.1 x 10~ 20 s. Zero friction (/ = 0) has been assumed. 



rather weak. 

We would like to point out that the theoretical values for Pl obtained in this 
work are smaller than those reported in our earlier calculations ||, |). There are 
two reasons for this: (a) the size of the lattice and (b) the lattice representation of 
the first derivative operator in the Dirac equation. Because of constraints in the 
amount of computer time available to us we utilized a smaller cubic lattice in our 
prior calculations || with N x x N y x N z = 29 3 lattice points. More recently, we 
were able to increase the size of the lattice substantially, in particular in fission (z-) 
direction (see above). In Fig. 2 of ref. |ll| we have demonstrated the convergence 
of our results for the muon attachment probability in terms of the lattice size and 
lattice spacing. Another reason for the difference between the current and prior 
results is the lattice representation of the first derivative operator, Eq. (|l0|), in the 
Dirac equation. In ref. ||,|5| we utilized a combination of forward and backward 
derivatives for the upper and lower spinor wave function components; after extensive 
testing of Coulomb potential model problems with known analytical solutions we 
have found that the symmetric derivative operator provides a more faithful lattice 
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Fig. 2. Contour plot of the logarithm of the muon probability density at t = 
1.1 x 10~ 20 s shows no evidence of muon ionization. 

representation. The results reported here and in ref. have been obtained 

utilizing the symmetric derivative prescription. 

5. Comparison of Theory with Experiment 

There are only a few experimental data available for comparison. Schroder et al. 
|jj"2f measured for the first time mean lifetimes of muons bound to fission fragments 
of several actinide nuclei. The muon decays from the K-shell of the muonic atom 
through various weak interaction processes at a characteristic rate A = Ao+A c , where 
Ao = (2.2 x 10 _6 s) _1 is the free leptonic decay rate for the decay process fi~ — > 
e~ + z7 e +j/ M and A c denotes the nuclear capture rate; A c depends upon the charge and 
mass of the fission fragment. From the observed lifetime r M = 1.30 x 10 _7 s Schroder 
et al. estimated an upper limit for the muon attachment probability Pl < 0.1. It 
must be emphasized that this number represents an integral over the whole fission 
mass distribution and, hence, cannot be directly compared to the numbers given in 
Fig. § 

The most complete experiments have been carried out by Risse et al. |l3| ] at the 
Paul Scherrer Institute (PSI) in Switzerland. The basic experimental approach is 
to place a fission chamber inside an electron spectrometer. The incident muons are 
detected by a scintillation counter. An event is defined by a (jjT , /i/2e — ) coincidence 
where the fission fragments are observed in prompt and the muon decay electrons 
in delayed coincidence with respect to the incident muon. The magnetic field of 
the electron spectrometer allows for a reconstruction of the electron trajectories. 
Thus, it is possible to determine whether the muon decay electrons originate from 
the heavy or the light fission fragment. 
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Fig. 3. Muon attachment probability to the light fission fragment as function of 
nuclear energy dissipation for 2 ||Np. Results are shown for fragment mass asym- 
metries £ = 1.05 (upper curve), 1.10, 1.15, and 1.20 (lower curve). 



For several mass bins of the light fission fragment, muon attachment probabili- 
ties Pi, have been measured; the experimental data are given in Table @. It should 
be emphasized that the mass bins are relatively broad. Because the theoretical 
values for Pl depend strongly on the mass asymmetry it is not justified to assume 
that Pl remains constant within each experimental mass bin. Instead, to allow for 
a comparison between theory and experiment, we have to multiply the theoretical 
Pl values in Fig. | with a weighting factor that accounts for the measured relative 
mass distribution |13l of the prompt fission events within this mass bin. We sub- 
sequently integrate the results over the sizes of the experimental mass bins. Due 
to the relatively low excitation energy in muon-induced fission, the fission mass 
distribution exhibits a maximum at £ = Ah/Al = 1.4 and falls off rather steeply 
for values larger or smaller than the maximum. This means that the large values of 
Pl ~ 0.5 at or near fission fragment symmetry £ = 1.0 will be strongly suppressed. 
The resulting theoretical values for Pl are given in the last column of Table |l| It 
is apparent that our theory agrees rather well with experiment. Because of the 
size of the error bars in the experiment and because of the weak dependence of the 
theoretical values of Pl on the dissipated energy, it is not possible to extract very 
precise information about the amount of energy dissipation. 

From a comparison of our theoretical result for the mass bin Al = 118.5 — > 
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Table 1. Muon- attachment probabilities to the light fission fragment, Pl, for 
237 Np(/i~, /). Exp. data are taken from ref. p3[ . 



mass bin Al 


mass asymmetry 


P L (exp) 


Pi(theo) 


118.5 -» 111.5 


1.000 -> 1.126 


(25.5 ± 8.5) X 10" 2 


26.0 x 10- 2 ,£ diss = OMeV 
22.3 x 10- 2 ,£ dias = 22MeV 


111.5 -* 104.5 


1.126 -> 1.268 


(9.7 ±2.6) x 10" 2 


6.62 x 10- 2 ,£ diss = OMeV 
3.51 x 10- 2 ,£ diss = 22MeV 



111.5 with the measured data we extract a dissipated energy of order 10 MeV for 
237 Np while the second mass bin A L = 111.5 — > 104.5 is more compatible with 
zero dissipation energy. We place a higher confidence on the theoretical results for 
the first mass bin because the probabilities Pj, are substantially larger and hence 
numerically more reliable. We like to point out that our theoretical value -Ediss = 10 
MeV is compatible with results from other low-energy fission measurements that 
are based on the odd-even effect in the charge yields of fission fragments p4|] . In 
addition to 237 Np we have also studied muon-induced fission of 238 U; the results for 
muon attachment are very similar p"l| ] . 

6. Conclusions 

We have studied the dynamics of a muon bound to a fissioning actinide nucleus 
by solving the time-dependent Dirac equation for the muonic spinor wavefunction; 
the fission dynamics is described classically. The theory predicts a strong mass 
asymmetry dependence of the muon attachment probability Pl to the light fission 
fragment; this feature is in agreement with experimental data. Our calculations 
show no evidence for muon ionization during fission. The theory also predicts a 
(relatively weak) dependence of P L on the dissipated energy. By comparing our 
theoretical results to the experimental data of ref. we extract a dissipated 

energy of about 10 MeV for 237 Np (see Table 1). Using the dissipation function 
defined in Eq. (||), this value corresponds to a fission time delay from saddle to 
scission of order 2 x 10 -21 s. 
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